UNIVERSAL BEHAVIOR OF THE COEFFICIENTS OF THE CONTINUOUS 
EQUATION IN COMPETITIVE GROWTH MODELS 
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The competitive growth models involving only one kind of particles (CGM), are a mixture of two 
processes one with probability p and the other with probability 1 ~ p. The p— dependance produce 
crossovers between two different regimes. We demonstrate that the coefficients of the continuous 
equation, describing their universality classes, are quadratic in p (or 1 — p). We show that the 
origin of such dependance is the existence of two different average time rates. Thus, the quadratic 
p— dependance is a universal behavior of all the CGM. We derive analytically the continuous equa- 
tions for two CGM, in 1 + f dimensions, from the microscopic rules using a regularization procedure. 
We propose generalized scalings that reproduce the scaling behavior in each regime. In order to 
verify the analytic results and the scalings, we perform numerical integrations of the derived analyt- 
ical equations. The results are in excellent agreement with those of the microscopic CGM presented 
here and with the proposed scalings. 
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Evolving growing interfaces or surfaces can be found 
in many physical, chemical and biological processes. For 
example, in film growth either by vapour deposition 
or chemical deposition P, 01, bacterial growth and 
propagation of forest fire The resulting interface 

has a rough surface that is characterized through scal- 
ing of the interfacial width W defined as W{L,t) — 
{{[hi — (^i)]^)^^^}, where hi is the height at the position 

i, (hi) — hi is the spatial average, L is the linear 

size, d is the spatial dimension and {} denote^configura- 
tion averages. The general scaling relation l] for these 
growing interfaces that evolves through a single model 
can be summarized in the form W(L,t) ^ L" f{t/L'^), 
where the scaling function f{u) behaves as f{u) ^ 
(/? = z/ct), for u <^ 1 and f{u) ^ const for u 1. The 
exponent a describes the asymptotic behavior where the 
width saturates due to finite sizes effects, while the expo- 
nent P represent the early time regime where finite-size 
effects are weak. The crossover time between the two 
regimes is tg = . 

The study of growth models involving one kind of par- 
ticles in competitive processes (CGM) has received little 
attention, in spite of the fact that they are more realistic 
describing the growing in real materials, where usually 
there exist a competition between different growing pro- 
cesses. As an example, in a colony of bacteria growing 
on a substrate, a new bacteria can born near to another 
and stay there, move into another place looking for food 
or died. This "bacteria" can be thought as a particle un- 
dergoing either a deposition/evaporation process or de- 
position/surface relaxation. 

The processes involved in the CGM could have differ- 
ent characteristic average time rate. Recently Shapir et 
al. 1^ reported experimental results of surface rough- 



ening during cyclical electrodeposition dissolution of sil- 
ver. Horowitz et al. introduced a competitive growth 
model between random deposition with surface relax- 
ation (RDSR) with probability p and random deposition 
(RD) with probability 1 - p, called RDSR/RD. The au- 
thors proposed that the scaling behavior is characteris- 
tic of an Edward Wilkinson (EW) equation, where the 
coefficient associated to the surface tension v depends 
on p. The dependance of on p governs the transition 
from RDRS to RD. Using a dynamic scaling ansatz for 
the interface width W they found that the results are 
consistent provide that v (x p^ . Also Pellegrini and Jul- 
lien 7] have introduced CGM between ballistic deposi- 
tion (BD) with probability p and RDSR with probability 
1— p, called BD/RDSR. For this model Chame and Aarao 
Ries |Sj presented a more carefully analysis in 1 -I- 1-d and 
showed that there exist a slow crossover from an EW to 
a Kardar-Parisi-Zhang (KPZ) for any p > 0. They also 
found that the parameter p is connected to the coefficient 
A of the nonlinear term of the KPZ equation by A ~ p''', 
with 7 = 2.1. 

In this letter, we show that the origin of such depen- 
dance is the existence of two different average time rates. 
Thus, the quadratic p— dependance is an universal fea- 
ture of all the CGM. To our knowledge this is the first 
time that the p-dependancc on the coefficient of the con- 
tinuous equations is obtained analytical from the micro- 
scopic dynamics. 

In order to test our hypothesis, we derive the an- 
alytical continuous equations for the local height for 
the RDSR/RD and BD/RDSR models. The proce- 
dure chosen here is based on regularization and coarse- 
graining of the discrete Langevin equations obtained 
from a Kramers-Moyal expansion of the master equation 
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Lets introduce first the general treatment of this prob- 
lem. Let us denote by hi{t) the height of the i-th 
generic site at time t. The set {hi, i = 1,...,L} de- 
fines the interface. Here we distinguish between two 
competitive processes: A with probability p and aver- 
age time of deposition ta, and B with probability 1-p 
and average time of deposition tb- In deposition pro- 
cesses with p = 1 the average rate of deposition is given 
by Tq"^ = {dhi / dt} ^p^j^y If the process is made with 
probability p the average rate of deposition is given by 
= {dhi / dt\ ^^-^ = p {dhi / dt} ^p^^y The same hold for 
a process with probability (1 — p). Thus, the particles 
are deposited at an average rate 
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In the average time of each process, the height in the site 
i increases by 



h^(t 



-TA) 

tb) 



h,{t) 



pRf, 
Ai-P) Rf 



(2) 



where Rf and Rf are the growing rules for processes A 
and B respectively and a_L is the vertical lattice spacing. 
Expanding hi{t -\- ta) and hi{t + tb) to second order in 
Taylor series around ta and tb , we obtain 



hi{t + Tj) - h,{t) 



dhj 
~dt 



Tj 



(3) 



for the process J — A,B. Thus, the evolution equation 
for the height (in the site i) for this CGM is given by 



dhi 
~dt 



(4) 



where the first moments of the transition rate for each 
process are 
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and the Gaussian thermal noise rji (t) has zero mean and 
covariance 

HitHit')} = + K^'-''^) % Sit - t') . (6) 

In order to test our analytical result, we use two mod- 
els. The first model RDSR/RD considers a mixture j^] 
of RDSR (process A) with probability p and RD (process 
B) with probability 1—p. Lets introduce the growth rule 
for each process for the first model. In the RD growth 
model one chose a column of a lattice, at random, among 
L and a particle is launched until it reaches the top of 
the selected column. The RDSR is a variant of the RD: 
a particle is released from a random position but when 
it reaches the top of the selected column is allowed to 



relax to the nearest neighbor (nn) column if their height 
are lower that the selected one. If the height of both 
of the nn are lower than the selected one the relaxation 
takes place with equal probability to one of them. For 
RD, W{L, t) does not depend on L, this means that the 
width W does not saturate due to the lack of lateral cor- 
relations. Thus, in this model: W{t) ~ Moreover, 
the RDSR model generates lateral correlations, therefore 
one has (inDSR = 1/4 and urdsr = 1/2. The first mo- 
ment of the transition rate for these processes are 



ft-^l'^) - ,(2) , , ,(3) (4) 
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(1) where the rules for both processes can we written as 
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Q{H^^) mr') , (8) 
{i [1 - mi^')\ + mi^')} [1 - mr')] , 
{i [1 - mr')\ + mr')} [i - ■ 



where H^^^ = {hi±s — hi±k)/a±, and Q{z) is the unit 
step function defined as 8(z) = 1 for z > and 8(z) = 
for z < 0. The representation of the step function can 
be expanded as Q{z) — J2T=o'^kz'' providing that z is 
smooth. In any discrete model there is in principle an in- 
finite number of nonlinearities, but at long wavelengths 
the higher order derivatives can be neglected using scal- 
ing arguments, since one expect affine interfaces over a 
long range of scales, and then one is usually concerned 
with the form of the relevant terms. Thus, keeping the 
expansion of the step function to first order in his argu- 
ment and replacing the expansion Eq. (|HJ| , Eq. Q can be 
written as 



dhi a±{l-p) aj_p f A'^hi 

— = + 1 + ci 

dt TB TA \ a_L 



Vrit) , (9) 



where A'^hi — hi+i — 2hi + hi^i ~ a| d'^h/dx^\hi, and 
ay is the horizontal lattice spacing. Replacing the rates 
given by Eq. in Eg^ @ and using a standard coarse- 
grain approach [T]| the continuous equation for this 
CGM is 

dh , ^ , ^ d'^h , ^ , ^ 

— =F{p)+v{p) — +ri{x,t), (10) 



where h = h{x,t) and 



Fip) = ^[il-pf+p^] 

To 

v{p) = 2ci^p2 

To 



(11) 



The noise covariance is given by 

Mx, t)r,ix',t')} = Dip) S{x - x')5{t - t') , (12) 
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where 

D{p) = ail F{p) . (13) 

Equations Ullfl and Eq. (|13|l shows that the quadratic de- 
pendance on the coefRcients of the continuous equation, 
arises naturally as a feature of the CGM and is due to 
the existence of different average time rates. 

The second model is a mixture of RDSR with prob- 
ability 1 — p {B process) and ballistic deposition (BD) 
with probability p {A process) Q. The evolution rules 
for RDSR are ojf, with j = 2, 3, 4 [see Eq. ©]. In the BD 
model, the incident particle follows a straight trajectory 
and sticks to the surface at the column i. The height in 
the column i is increased in max[/ii + I, hi^i, hi-i]. If 
this process is done with probability p (A process), the 
rules can be summarized as: 

^ = 0(^+i) 0(^-i) , (14) 

+ Hi-' [i-mu)]} , 

where S{z, 0) — 0(z) + Q{~z) — 1 is the Kronecker delta. 
Following the steps leading to Eq. (|10(l the evolution 
equation for this process can be written as: 

^.FM + .(ri0 + A(ri(|)%„..O ,15) 

where 

F{p) = ^[il-pf + 4p^] , 

To 

v{p) = -i [ip2(i_co-2coCi) + 2ci(l-p)2] ,(16) 

Kp) = — —P ci (5 - 4co - ci) . 

To a±_ 

The covariance of noise and D{p) is given by Eq. (|12|l and 
Eq. H13|) , respectively. Notice that we have to change p by 
1 — p in all the above equations for RDSR, because in the 
first model RDSR is a kind A process and now is a kind 
B process. Equation (|16|l shows again that quadratic dc- 
pendance on the coefRcients of the continuous equation. 
The quadratic dependance of A on p, found by Chame 
and Aarao Reis is a general feature of the CGM. 

As both models have an EW behavior, it is expected 
that in that regime the following generalized scaling 
ansatz lEH, 

W\p, L, t) ^ L^" [D{p)/i.{p)] f {iy{p) tjU-) , (17) 

where f{u) ~ u^^ for m ^ 1 and f{u) ^ const for u ^ 1. 
Moreover, the second model is represented by a mixture 



of EW and KPZ universality classes. In the early time 
regime W{t) ^ i^rdsr ^ while a crossover to a KPZ, with 
Pkpz = 1/3 and a = axpz — 1/2, is expected in the 
intermediate regime before the saturation. Thus, for the 
KPZ regime we propose the following generalization 
of the scaling behavior of the width 

W^{p,L,t)^ L^'^[D{p)/v{p)] f(x{p)^D{p)/i,{p)t/L^) , 

(18) 

where z = 3/2, and f{u) ~ u'^Pkpz for u < 1 and f{u) ~ 
const for w ^ 1. 

In order to test our analytical result and the proposed 
scalings, we perform a numerical integration of Eq. H10(l 
and Eq. \ I ^lotice 

lo'l l-r™ 1-™ 1 
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FIG. 1: (a) Log-log plot of W'^v(j>) / D(p) for the RDSR/RD 
model as function of v(p) t for L — 128. The different symbols 
represent different values of p, p = 0.04 (o), p — 0.08 (□), 
p = 0.016 (A), p = 0.32 (+), and p = 0.64 (*). Here we 
used C — 2.58 and b — 0.2 as parameters of the 0-function 
representation. The dashed lines are used as guides to show 
the RD regime with 2/3 = 1 and the EW regime with 2/3 = 0.5. 
(b) Log-log plot of W^iy{p)/D{p) for the BD/RDSR model as 
function of v{p) t for L — 1024. The symbols represent the 
same as in Fig. (la). Here we used C — 0.18 and b — 0.5. The 
collapse of the curves at the earlier stage clearly shows the 
EW behavior (2/3 — 0.5). After this stage the curves split and 
undergoes a slow crossover to the KPZ behavior (2/3 — 0.66). 
The dashed lines are used as guides to show the EW regime 
with 2/3 — 1/2 and the KPZ regime. The slope showed here 
is 2/3 = 0.61. 

that in order to numerically integrate the continuous 
equation, we do need a continuum representation of the 
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FIG. 2: Log-log plot of as function of i as defined in the 
text, for p = 0.16 (o), p = 0.32 (□) and p = 0.64 (A). The 
empty symbols correspond to L = 512 and the filled ones 
to L = 1024 . The collapse of the curves on the saturation 
regime using z = 3/2 shows that the curves saturates with a 
KPZ behavior as expected. 

G-function to numerically compute the coefficients cq and 
ci related to the ones of the continuous equations. To 
perform the numerical integration, we chose the shifted 
hyperbolic tangent 16] as the continuous representation 
of e-function defined as Q{x) = {1 + tanh[C(a; + b)]}/2, 
where b is the shift and C is a parameter that allows to 
recover the Q in the limit C — > oo. The numerical in- 
tegration was made in short lattices using a discretized 
version of the continuous equations Eq.l^UJ and Eq. (fT^ . 
The results in large systems and the details of the inte- 
gration are beyond the scope of this letter and will be 
published elsewhere. 

For the first model, Horowitz et al. presented 
their data from simulations plotting the scaling rela- 
tion W/L^p-^ vs t/L^'p-y. Clearly, their (5 = 1 and 
y — 2 is related to our i/{p) and D{p) [see Eq. (|ll|l and 
Eq. (O]. In Fig. Ula) we plot W'^ iy{p) /D{p) as func- 
tion of v{p) t for the that model for different values of p 



and L = 128. This figure represent the same as in 
after coarse-graining. The agreement with the results of 
our numerical integration, the numerical simulation 
and the scaling presented in Eq. p7l) is excellent. On 
the other hand, for the second model, Chame and Aarao 
Reis [si did not present the result for W. They studied 
the crossover from EW to KPZ using an indirect method 
because of the slow convergence of the discrete model to 
KPZ behavior. The crossover is well represented in our 
Fig. nib) , where we plot the same as in Fig. QJa) but for 
the second model. It is clear the collapse of the curves in 
the EW regime. In the intermediate regime the KPZ be- 
havior appears thus, it is expected that Eq. (|18() holds in 
that regime. In Fig.Hwe plot W'^ = v{p)/[L^°'D{p)] 
as function of i = \{p)yj D{p) /v{p) t/L^ for three differ- 
ent values of p using ^ = 3/2. As p increases, the KPZ 
behavior appears earlier, but independent of p all the 
curves saturate as a KPZ. The agreement with Eq. H18|) 
is excellent in the saturation regime. The departure in 
the intermediate regime is due to a slow crossover to the 
KPZ and to finite size effects. 

Finally, notice that the quadratic dependence of the 
coefficients of the continuous equation on p is indepen- 
dent of the CGM considered, because it is due to two 
different rates of deposition given by Eq. This de- 
pendence is totally generally, as shown from Eq. |^ to 
Eq. 0. 

In summary, we demonstrate that the coefficient of 
the continuous equation have quadratic dependance on p. 
This feature is universal for all the CGM model and is due 
to the competition between different average time rate. 
We propose generalized scaling for the model that repro- 
duce the scaling behavior in each regime. The numerical 
integration of the continuous equation are in excellent 
agreement with the propose scalings and the numerical 
simulation of the models. 
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